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ABSTRACT 

The return to the quiescent state of the Anomalous X-ray pulsar XTE J1810— 197 following 
its 2003 outburst represents a unique opportunity to probe the surface emission properties of 
a magnetar The quiescent emission of XTE J1810— 197 is composed of two thermal compo- 
nents, one arising from the whole star surface, and the other from a small warm spot on it. 
By modeling the magnitude and shape of the pulse profile in narrow spectral bands, we have 
been able to constrain the physical characteristics and geometrical parameters of the system: 
the two angles that the line of sight and the spin axis make with respect to the warm spot 
axis {ifj and ^ respectively), the angular size of the spot, and the overall surface temperature 
distribution. Our modeling accounts for the general relativistic effects of gravitational redshift 
and light bending near the stellar surface, and allows for local anisotropic emission. We found 
that the surface temperature distribution on the neutron star is consistent with the expectations 
of a dipole magnetic field configuration; the local radiation requires a pencil-beamed emis- 
sion pattern, suggesting the presence of a magnetized atmosphere. For a typical value of the 
radius, R=13 km, the viewing parameters (symmetric for an interchange between ij} and ^), 
range from i^ — ^ = 38° to {ip, ^)=(52°,29°). These angles are consistent with those obtained 
by modeling the AXP in outburst, with uncertainty contours reduced by a factor of 2.5. 

Key words: pulsars: general - stars: neutron - X-rays: individual: XTE J1810— 197. 



1 INTRODUCTION 

XTE J1810— 197 is an isolated neutron star (NS) belonging to the 
class of the Anomalous X-ray Pulsars (AXPs); these objects, to- 
gether with the Soft Gamma-ray repeaters (SGRs), are believed to 
be magnetars: isolated neutron stars whose thermal emission and 
occasional outbursts are powered by their extremely strong mag- 
netic fields (Duncan & Thompson 1992; Thompson & Duncan 
1995). Originally one of the thousands of the faint X-ray sources 
cataloged by ROSAT, the variable nature of XTE J1810-197 was 
revealed by the outburst of 2003, with its sudden increase in X- 
ray luminosity by a factor ~ 100, decaying on a time-scale of 
years. The source rotation period and its derivative were conse- 
quently measured and found to be P = 5.54 s and P = 1.1 — 
2.1 X 10~^^ s/s. These timing properties imply a magnetic field 
Sdip ~ 3 X 10^^ G, confirming the magnetar classification of 
the source (Ibrahim et al. 2004, Gotthelf at. al. 2004). The source 
was monitored repeatedly for more than 7 years with several X-ray 
observatories, up to the return to quiescence (Gotthelf & Halpern 
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2005, 2007; Halpern & Gotthelf 2005; Bemardini et al. 2009; Al- 
bano et al 2010). These studies presented a unique opportunity 
to probe the emission mechanisms of a strongly magnetized NS 
by taking advantage of the flux evolution during its decay. While 
analysis of phase-averaged spectra alone cannot uniquely distin- 
guish among competing emission models, the addition of the steady 
change of the spectrum and pulse profile over time greatly increases 
the diagnostic power. 

Perna & Gotthelf (2008) developed a detailed emission model 
for the energy-dependent pulse profile of XTE J1810— 197 follow- 
ing its outburst. This model, which was tailored to the specific 
surface emission distribution in the post-outburst phase, can take 
on any viewing geometry, includes the general relativistic effects 
of light deflection and gravitational redshift, and allows for local 
anisotropic emission. The application of this model to the first 4 
sets of XMM-Newton data acquired during the temporal evolution 
of the flux from XTE J1810— 197 following the outburst (Septem- 
ber 2003 - September 2004), provided a constraint on the under- 
lying emission geometry and radiation properties of this transient 
magnetar in its post-outburst phase. 

In this paper we present the results of the modeling of the 
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Figure 1. XMM (red) and Chandra (blue) 0.5-10 keV flux measurements 
of XTE J1810— 197 following its outburst. It is evident from the data 
that XTE J1810— 197 has now reached quiescence. The dashed line rep- 
resents the X-ray flux level (~ 7.5 X 10"^'^ ergcm^^ s^^) as recorded 
by ROSAT, Einstein, and ASCA before the outburst onset (Ibrahim et al. 
2004, Gotthelf et al. 2004). The analysis presented in the paper is performed 
over three close XMM-Newton pointings represented by the last data point. 



spectral and timing data of XTE J1810— 197 obtained with 3 com- 
bined XMM-Newton pointings (September 2009) upon the return 
of the source to quiescence (see Figure [T] last point). Indeed, our 
goal here is that of studying the properties of the quiescent emission 
of this magnetar, which carries information on the surface tempera- 
ture distribution of the star, and hence on its magnetic topology. We 
model these data using a modified version of the emission model 
by Pema & Gotthelf (2008), updated to include: (a) the presence of 
quiescent emission from the full surface of the star, (b) the changed 
and reduced emission from the region heated by the outburst. 

The new data and its spectral and timing analysis is reported 
in §2. In §3, we discuss the properties of the quiescent emission. 
The theoretical emission model is described in detail in §4, and the 
results of its application to the data are given in §5, following with 
a discussion in §6. 



2 OBSERVATIONS 

The quiescent state of the source was observed by XMM-Newton 
during three consecutive close pointings (18 days total time span) 
in 2009, on September 5, 7, 23 for 19, 18, 12 ks, respectively. All 
the observations were performed with the PN instrument (Striider 
et al. 2001) in large window mode, with the medium filter applied, 
and the MOS 1 and MOS 2 instruments (Turner et al. 2001) in small 
window mode, with the use of the medium and thin filter, respec- 
tively. The PN resolution time with this configuration is 47.6 ms 
while the MOS 1/2 resolution time is 0.3 s. Data were processed 
with SAS version 10.0.0, using the updated calibration files (CCF) 
available in August 2010. Standard data screening criteria were ap- 
plied in the extraction of scientific products. Time window crite- 
ria were used for removing time intervals contaminated by solar 
flares. A total of 48 ks of good exposure time was obtained. Pho- 



ton arrival times were converted into barycentric dynamical times 
(TBD) using the SAS tool barycen and the milliarcsec radio po- 
sition of Helfand et al. (2007): RA=18'' 09" 51=0870, Dec=-19° 
43' 5l".931 (J2000). 

Source photons were extracted from a centered circular region 
of radius 55" containing 90% of the source counts. For both the 
timing and spectral analysis we extracted the background from the 
same PN or MOS CCD where the source lies using a circular re- 
gion of the same size as that of the source. In the following analysis, 
we combine the three PN spectra collected by XMM-Newton over 
an interval of 18 days, after verifying that the individual spectra 
were consistent with each other within the uncertainties. The spec- 
tra were binned to have at least 30 counts per fitting channel to 
insure adequate fit statistics. 

For each observation and instrument we extracted and 
summed data into four energy bands (following Pema & Got- 
thelf, 2008): 0.5-1, 1-1.5, 1.5-2 and 2-3 keV. Lightcurves were 
generated at the minimum allowed bin time of 0.3 s. We phase 
connected the three multi-instrument lightcurves using the phase- 
fitting technique outlined in Dall'Osso et al. (2003). The timing 
solution, referred to epoch MJD 54079, includes only one term, 
the rotation period, which was found to be P = 5.5406556 ± 
2 X 10^^ s. The brief span of the observations did not allow for 
a spin-down measurement, but only for a 3cr upper limit of P ^ 
9 X 10"^^ s/s (consistent with the value P = 0.8 - 1.0 x 10~" s/s 
reported by Camilo et al. 2007 for the 2006-2007 observation pe- 
riod). The pulse profile at different energy bands was determined 
by means of the P value reported above. The profile was found to 
be nearly sinusoidal and energy independent in analogy with earlier 
epochs (Halpem and Gotthelf 2005; Bernardini et al. 2009). 



3 PROPERTIES OF THE QUIESCENT EMISSION 

Halpern & Gotthelf (2005) and Gotthelf & Halpem (2005) showed 
that the post-outburst spectrum was composed of a multi-blackbody 
(BB) made up of two thermal components, which they interpreted 
as a warm ring surrounding a hot spot. Bernardini et al. (2009) fur- 
ther showed evidence for a third, cooler thermal component consis- 
tent with emission from the whole surface of the star. The recorded 
outburst flux was a factor ~ 100 times higher than in the quiescent 
state. They also found that, while the intensity of the emission from 
the whole NS surface was constant during the outburst (and equal 
to the one recorded in quiescence), the flux from the two hotter re- 
gions decreased exponentially with time (on a timescale r ~ 1 yr). 
Spectral analysis showed that the warm and the hot regions were 
shrinking with time, and the total luminosity of the star was con- 
sequently declining towards the quiescent level. When the source 
eventually returned to quiescence, its spectrum was fully consistent 
with that recorded by ROSAT before the outburst. More in detail, 
the quiescent flux was found to be composed of a cool BB compo- 
nent consistent with emission from the whole NS surface (for the 
best estimated distance value of ~ 3.3 kpc; Camilo et al. 2006, 
Minter et al. 2007, Durant & van Kerkwijk 2006), and a warmer 
BB emission coming from a small, residual spot. The quiescent 
spectrum of the source is therefore composed by two BB only (see 
Figure[2]left panel). It should however be noted that the quiescent 
spectrum also displays a statistically significant (~ 5cr) absorption 
feature around 1.1 keV. This feature, discussed in Bernardini et al. 
(2009), is of unknown origin. It could be a proton cyclotron line 
if the magnetic field is P = 2.2 x 10^* G (an electron cyclotron 
line would imply a B field about 2000 times weaker, out of the 
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Figure 2. Left panel: The quiescent spectrum of XTE 11810—197 modeled with a 2BB model with the softer component associated with the entire surface 
and the hotter component, with a higher level of PF, associated with a localized warm spot. Right panel: Pulse profile of the surface component (black square) 
compared to the warm spot pulse profile. Both profiles show a peak at the same phase interval (cji ~ 0.85). 



magnetar range). In our spectral fits, the feature is modeled with 
an edge. We have verified that, within uncertainties, the absorption 
depth {tc) at the energy threshold (Ec) is independent of phase, and 
of magnitude < Tc >= 0.32 ± 0.02. Hence the presence of this 
feature is not expected to affect the results from the timing analysis 
and is not included therein. Such a feature is however included in 
all spectral fits in the present work. 

Since XTE J1810— 197 has now returned to quiescence, the 
temperature distribution on the star surface reflects the overall mag- 
netic field distribution, as the conductivity is enhanced along mag- 
netic field lines. Eor most AXPs, the high level of pulsed fraction 
cannot be produced by a temperature distribution following a dipole 
magnetic field (De Deo, Psaltis & Narayan 2000). In the case of 
XTE J1810— 197, on the other hand, the low level of pulsed frac- 
tion of the soft X-ray bands could be the result of such a field con- 
figuration, which we hence adopted as our starting point. A deter- 
mination of the XTE J1810— 197 surface temperature distribution 
would allow to unveil its magnetic field configuration. 

The spectral component corresponding to the cold NS surface 
easily dominates the emission below 1 keV, while the emission cor- 
responding to the warm spot dominates over 2 keV (see figure |2] 
left panel). Our first goal was the localization on the star surface of 
the warm spot with respect to the maximum intensity of the cooler 
surface emission (corresponding to the magnetic axis). Note that, 
albeit the NS surface is expected to have a temperature gradient, 
however, due to the limited S/N ratio of the spectral data, it can only 
be identified as a single BB in the spectral analysis (i.e. the colder 
one). The pulse profiles of the two BB components were conse- 
quently generated separately by taking into account only photons 
emitted in the 0.5—1 keV and 2—3 keV energy bands respectively. 
The result of this test shows (see figure[2]right panel) that the max- 
imum of the pulsed emission coming from the NS surface remains 
in phase with the maximum in the hardest energy band, dominated 



by the warm spot. Therefore, one the two (warmer) regions on the 
star surface associated with the magnetic poles must be very close 
to the center of the warm spot. It is important to note that, while the 
highest energy bands (dominated by the emission from the single 
hot spot), naturally produce produce a single-peaked profile, in the 
lowest energy band, dominated by the surface emission (character- 
ized by two symmetrically opposed warmer regions), a single peak 
profile is expected only for geometries for which only one pole is 
visible as the pulsar rotates (e.g. Page 1995). 

These results and considerations represent the basis to develop 
a model that allows us to predict the energy-dependent pulsed frac- 
tion and pulse profile, and use it to determine, for different values 
of the NS radius, the viewing geometry and beaming pattern of the 
emitted radiation that best match the observations. The model is 
described in detail in the next section. 



4 MODELING THE SPECTRUM AND THE PULSE 
PROFILE OF XTE JI8I0-I97 IN QUIESCENCE 

Following the motivations given in the previous section, the local 
temperature distribution, rth(r, ^) (in spherical coordinates), on 
the surface of XTE J1810— 197 is modeled as expected from ther- 
mal cooling under the influence of a dipole magnetic field (Heyl & 
Hernquist 1998; see also Perna et et al. 2001): 



Tth(6',0)=Tj, 



4 cos e^ 



3 cos^ 9p + 1 



(0.75 cos^ eip + 0.25)°-^ 



(1) 



where Tp is the pole temperature and Op is the angle between the 
radial direction at position [6, (j)) on the surface of the star and the 
magnetic pole. 

Superimposed to the emission from the whole (cold) surface 
of the star is the emission from a small (hot) spot, whose axis coin- 
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cides with the dipole axis. The mathematical description of the spot 
on the rotating surface (of temperature Th and angular radius Ph) 
follows the theory developed by Pechenick et al. (1993) with some 
generalizations presented by Perna & Gotthelf (2008). If ^ defines 
the angle between the spot/dipole axis and the rotation axis, and 
tp the angle between the observer's direction and the rotation axis, 
then the angle a that the axis of the hot spot makes with respect to 
the line of sight to the observer can be written as 



a(t) — arccos(cos-i/;cos^ + sirn/)sin^cos7(i)) 



(2) 



This angle is a function of the phase angle ■y{t) = Q{t)t swept by 
the star as it rotates with angular velocity Q{t). The surface of the 
star is described by the angular spherical coordinates {6,4>), and 
the coordinate system is chosen so that the z axis coincides with 
the direction of the line of sight to the observer. The hot spot is 
described by the conditions: 



^/3h 



if Q = 



and 



- /3h < 6» s; a + /3,, 



27r 



where 



€</}€ 



if a 7^ and Ph ^ a 



cos /3h — cos a cos 8 



sin a sin 6 
On the other hand, it is identified through the condition 



^e.(a,/3h, 



if Q 7^ and /3h > a 



(3) 



(4) 



(5) 



(6) 



where the outer boundary 6!^{a, j3h, (t>) of the spot is computed by 
numerical solution of the equation 

cos Ph ~ sin Of sin a cos (ji + cos 6*, cos a . (7) 

Due to the strong NS gravitational field, photons emitted at the NS 
surface suffer substantial deflection on their way to the observer. 
A photon emitted at a colatitude 8 on the star makes an angle 5 
with the normal to the surface at the point of emission. The relation 
between 5 and 6 is given by the ray-tracing functior|j (Page 1995) 



6(5) = 



Rs/2R 



X du 



-^ 



R 
2R 



^ _(l_2uKx2 



having defined x = sin 5. Here, R/Rs is the ratio between the NS 
and the Schwarzschild radius, Rs — 2GM/c? (we assume M = 
1.4Mq). 

A blackbody model for the local emission is assumed. While 
it would be desirable to perform this analysis with realistic mag- 
netized atmosphere models, the lack of an extensive set of such 
models for high B-field strengths, {B ^ 10^* G) and arbitrary 
inclinations (with respect to the NS surface) makes this more com- 
plete analysis not yet possible. This is particularly the case for the 
present study, since we are modeling the emission from the en- 
tire surface of the star, and over this there are large regions with a 
non-normal B. We note, however, recent work extending NS emis- 
sion models to non-normal fields. In particular, LLoyd (2003a,b) 
presented model spectra for B ^ 10^^ G and for arbitrary orienta- 
tion, for pure Hydrogen composition, and in the limit of complete 
ionization. Ho, Potekhin & Chabrier (2008) constructed partially 
ionized Hydrogen models for arbitrary field orientation and for 

^ In the emission code, to improve the computational efficiency of the 
above equation, we use the approximation derived by Beloborodov (2002). 



strengths in the range lO" < B ^ 3 x 10^^ G. While these models 
are very useful for exploring NSs with moderate fields (e.g. Mori 
& Ho 2007), they are still not appropriate for the magnetic field 
strengths needed for a self-consistent modeling of XTE 11810-197 
(Bdip of about 3 x 10^* G). Hence here we adopt the empirical ap- 
proach of parameterizing the level of anisotropy with the function 
f{5) oc cos" 5 (since 'pencil' beaming dominates in atmosphere 
models), and perform the timing analysis with different values of 
n within a reasonable range as suggested by beaming in realistic 
descriptions of magnetized atmosphereo 

The observed spectrum as a function of phase angle 7 is then 
obtained by integrating the local emission over the observable sur- 
face of the star, including the effect of gravitational redshift of the 
radiation (Page 1995) 



F{E^,-l) 



27^RL 



^2g-iVH.(B^) 



2xdx 



^Ioie,^)n[E^e-^^-Tie,<t>)] 



(9) 



in units of photons cm^^ s~^ keV~^ . In the above equation, the ra- 
dius and energy as observed at infinity are given by R^c ~ Re^ " , 
and Eao = Ee ' , where R is the star radius, E is the energy emit- 
ted at the star surface, and As is defined through the relation 



R 



(10) 



For the spectral function, given by n{E, T) — 
l/[exp{E/kT) - 1], the temperature T{e,(j)) is equal to Th 
if {8, (f)} satisfy any of the conditions l[3j through (O, and it is 
given by T{8,(j)) = Tth{8,(l)) otherwise. Correspondingly, the 
weighted intensity Io{d,(j)) is given by the beaming function 

fimi 

The phase- averaged flux is then readily computed as 



1 f^'' 

Favol-Boo) = — / d'yF{E^,'y) 
^1" Jo 



(11) 



Note that the phase dependence 7 in Eq.{9]l comes from the viewing 
(8) angles implicit in a{t) and from the series of conditions l|3j through 
Q. 

We also included a multiplicative factor which accounts for 
the hydrogen column density between the observer and the star. 
Note that absorption does influence the predicted pulsed fractions, 
when these are computed over finite energy intervals (Perna et al. 
2000). The magnitude of absorption was fixed to that obtained from 
the spectral fit (see i]5.1t . The model described in this section was 
imported into the spectral fitting package XSPEC (Arnaud 1996), 
and used to fit both spectra and pulse profiles. 



^ In the magnetized atmosphere models by van Adelsberg & Lai (2006), 
the beaming strength depends on the magnitude of the B field, on the at- 
mosphere temperature, on the observation energy, as well as on the angle 
5 itself. For fields and temperatures in the magnetar range, the models by 
van Adelsberg & Lai (2006) predict a strong forward beaming for angles 
"5 S; 40 — 60 deg, and a much lower anisotropy level (fan-like) at larger 
angles. For example, for B ~ 10^** G and T ~ 0.4 keV, an approximation 
to the intensity for 5 < 50 deg is f{S) oc cos*^-* 5 at i? ~ 0.2 keV, and 
/(<5) oc cos^ 5 at £ ~ 1 keV. 
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5 RESULTS 

5.1 Spectral modeling 

Fitting the phase averaged-spectrum with the model in Eg. Jilt re- 
quires fixing the viewing angles ^,1/'. However, since the viewing 
geometry is not known a priori, we followed an iterative proce- 
dure for our combined spectral and timing analysis. We start by as- 
suming 7/) = 5 = 0, the face-on geometry, with the observer look- 
ing directly down the co-aligned rotation axis and magnetic pole. 
This implies a = and speeds up the spectral fitting substantially 
as the geometric condition describing the hot spot (cfr. Eq. [Sj de- 
pends only on the coordinate 6, reducing the flux integral in Eq.(|9j 
to a single dimension. After fitting the phase-averaged spectrum, 
the spectral parameters are used to compute the pulse profiles for 
5, ■(/>; this restricts the range of viewing angles ?/) and ^ to those that 
best match the pulsed fractions. We then refit the phase-averaged 
spectrum to refine the spectral parameters, but now use the best 
ip, ^ from the timing analysis. This procedure is iterated until the 
change in the parameters are consistent with the measurement er- 
rors. This procedure is found to converge in only one iteration and 
thus allows us to explore a wide grid of viewing angles (like e.g. in 
Gotthelf et al. 2010) in a reasonable time. 

To begin with, the source distance was first fixed to D=3.3 kpc, 
based on radio pulse dispersion measure (Camilo et al. 2006), and 
consistent with the measurement derived from HI absorption (D = 
3.5 ± 0.5, Minter 2007) and the measurement derived from Red 
Clump Stars in the direction of the source (D = 3.1 ± 0.5 Durant 

6 van Kerkwijk 2006). For completeness, we then also studied the 
cases corresponding to the distance fixed at 2 kpc and 5 kpc, which 
correspond to the borderline values of the 3a uncertainty on the 
distance. All the fits were performed in the 0.5—3 keV spectral 
band. Above 3 keV, source detection is not significant. We repeated 
the spectral analysis for several values of the beaming parameter n 
(in the range — 2.5, using steps of 0.5), and found that, within the 
measurement uncertainties, the inferred spectral parameters were 
all consistent. Hence, for the spectral analysis, we set n = 0. 

Like in the case of the analysis performed by Pema & Gotthelf 
(2008), an important technical issue for these fits is the degree of 
degeneracy between the radius R and the variable parameters, in 
this case (kTh, Ph, kTp). Without fixing the radius and the distance 
there is no unique solution, and we considered a range of possible 
values, 9 ^ -R ^ 15 km, in 1 km increments, for the radius, and 
the above 3a range of 2 ^ _D ^ 5 kpc. 

For the distance values of _D = 3.3 kpc and D — 2 kpc, over 
the sampled range of radii, spectral analysis provides acceptable 
fits, but did not allow a preferred radius based on the x^ measure- 
ment (see table[Tll. In the case of D = 5 kpc, spectral fits showed a 
significantly higher x^ value, ranging from 1.4 (for 7? = 15 km) to 
2.04 (for R = 9 km). These fits resulted to be statistically unaccept- 
able; consequently, the timing analysis described in the following 
is only performed for D = 3.3 kpc and D = 2 kpc. 



5.2 Pulsed fraction and pulse profile modeling 

Given the smooth and nearly sinusoidal pulse shape, the pulsed 
fraction (PF) of the signal has been determined using the expres- 
sion 



PF: 



F ■ 

-t niin 



(12) 



Fluxes are integrated over the given energy bands. In the adopted 
geometry, the maximum and minimum fluxes of the model, -Fmax 



and -Fmin, correspond to phases 7 = and j — tv, respectively. As 
discussed in §3, this is true for any combination of viewing angles 
in the higher energy bands ^ IkeV, in which the emission is either 
dominated, or largely influenced, by the (single) hot spot. In the 
lowest energy band (0.5- Ike V), where the contribution from the 
surface emission is dominant, a single peak can only be obtained for 
viewing angles ^,ip ^ (,,nax,^max- The latter depend mildly on 
n, R, and are generally < 50° —55°. Our viewing parameter search 
is restricted to the range of angles for which the lowest energy band 
remains single peaked. Incidentally, the high ^, tp parameter range 
is also independently ruled out by the PF in the highest energy band 
alone: for most combinations of n, R, the predicted PF would be 
much higher than the measured value. 

A measure of the PF in each energy interval was obtained us- 
ing the timing solution, phase-connecting the three pointings. The 
background level, which was variable during the three pointings, 
was subtracted for each multi-instrument lightcurve. We found 
that the PF increases with energy: PFjo.s-ikoV) = 17 ± 1%, 
PF(i_i.5kcV) = 26 ± 0.8%, PF(i.5-2kcV) = 36 ± 1.3%, 
PF{2-3kcV) = 47 ± 2.6% (all the uncertainties are hereafter re- 
ported at 1(T confidence level, c.L, unless otherwise stated). Above 
3 keV, due the low S/N ratio, only an upper limit on the PF could 
be obtained (but this is not an useful value for the purpose of this 
work). Hence our timing analysis uses the data up to 3 keV. 

Starting with the best spectral fit model parameters presented 
in Table [T] obtained for 1/; = ,^ = 0, the best values of ijj* and ^* 
needed to reproduce the observed pulsed fraction were searched (in 
the case of D=3.3 kpc and D=2 kpc). Given the high S/N of the data, 
both the pulsed fraction and the full pulse profile were modeled. 

For each value of the NS radius fitted for in i]5. 11 (9-15 km), 
the pulsed fraction, defined in Ea.il2t. was computed over the grid 
of angles {i,ip) < (Cmax.V'max) dcg, lu 1 degree intervals. For 
each value on the grid, the model predictions were compared with 
data. Notice that the flux depends on the angles ^ and ip only 
through the parameter a in Eq.Q, and therefore it is symmetric 
with respect to an exchange of ^ and i/;. 

In order to explore the behaviour of the PF with both R and n 
independently, we first examined how the PF varied with R at fixed 
71, and then how they varied with n at fixed R. In the following, we 
first report the results for n = 1 and i? = 9 — 15 km, in steps of 
1 km and then for i? = 9 km and i? = 15 km, for n = — 2.5, 
in steps of 0.5. The behaviour for other combinations of parameters 
can then be inferred from the results shown. 

For the case n — 1, Table |2] reports the best geometry (ex- 
pressed through the angles tp* and ^*), for each sampled radius, 
and the corresponding predicted PF as a function of the energy 
band, in comparison with the data. For this case, we further ex- 
plored the dependence of the results on the source distance. In the 
case of D = 3.3 kpc, low radius values (9 ^ -R ^ 12 km), un- 
derpredict the observed PF for low energy intervals 0.5 ^ -E ^ 2 
keV. Moreover, for R = 11 and R = 12 km, the model overpre- 
dicts the observed PF in the highest energy band (2 ^ _E ^ 3 keV; 
see[3]l. The different trend of the PFs with radius in the low and 
high energy bands is due to the fact that the flux in these is domi- 
nated by different components: the cooler component from the all 
NS surface in the low energy band, and the hotter small spot in the 
high energy band. The smooth temperature gradient of the former 
makes the PFs less sensitive to changes in the viewing geometry, 
and hence the dominant factor in determining the change in the PF 
with R is the general-relativistic suppression of the pulsed fraction 
as R decreases. On the other hand, the pulsed fraction produced by 
the small hot spot is much more sensitive to changes in viewing 
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Table 1. Spectral fit parameters as a function of the radius R. la c.l. uncertainty is reported. 
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Figure 3. Model predicted PF (for the best emission geometry as reported in Table|2j for different radius values, compared with data. Left panel is for D = 3.3 
kpc, while right panel is for D = 2.0 kpc. Beaming factor is equal to 1. Dashed and dot-dashed lines represent the model prediction for low radius values 
(9 ^ -R ^ 12 km), while continuous black lines are obtained for high radius values (13 ^ if ^ 15 km), la c.l. uncertainty is reported. 



angles with the result that smaller R's, which are best fit by larger 
values of ^, ip (see Table |2j, yield higher PFs. We found that it is 
not possible to reproduce the observed pulsed fraction for low ra- 
dius values (9 ^ i? ^ 12 km), which are hence rejected by the 
model for this distance. The predicted PF for high radius values 
(13 ^ R ^ 15 km) is instead fully consistent with the data, within 
Icr uncertainty. Figure [3] shows the model predicted PF, for all the 
values of the radii considered, compared with the data. Similar re- 
sults are obtained for D = 2 kpc; however in this case the range 
of allowed radii is wider, including also R — 12 km. For, e.g., 
an intermediate value of the radius, R = 13 km (D — 3.3 kpc), 
the best viewing geometry angles range from tp* = (^* — 38° to 
('tp,0 = (52° 29°) at 3cr c.l. (see Tabled. The angles {^ , that 
provide the best match to the PF data clearly vary with R, for a 



fixed value of the beaming factor n. Smaller radii require a larger 
variation in a{t) {a varies between tp — S, and ■i/; + ^) to compensate 
for the stronger general-relativistic suppression of the flux modula- 
tion. 

The xi map, computed for both the case of i? = 9 km and 
7? = 14 km (Figure Ell, displays the 68%, 90% and 99% confi- 
dence levels. This map is produced by comparing the model and 
observed pulsed fraction over a range of possible (^, tp) angle pairs, 
for our best fit spectral model parameters. The viewing geometry 
was found to be well constrained, with the range of allowed Sa so- 
lutions a factor of 2.5 smaller than (but consistent with) what was 
found by Perna & Gotthelf (2008) in the analysis of the outburst 
decay (where the 3a c.l. of the viewing parameters were ranging 
from tp* — ^* = 37° to tp,^ — 85°, 15°). This consistency is 
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Figure 4. Reduced chi-square (Xv) maps obtained by comparing data and 
model pulsed fraction described in the text for a range of viewing angles 
^ and g. The 68%, 90% and 99% confidence levels ai'e shown for the best 
match to the observed pulsed fractions using the beaming pattern n = 1 for 
R = 9 km and ij = 14 km. The results are clearly degenerate with respect 
to an interchange of (^, i/i). A comparison with Fig. 2 in Perna & Gotthelf 
(2008) shows the substantial reduction in the confidence range. 




Figure S. Predicted pulse profile for _R = 13 km (black line), and R = 
9 km (red line), D = 3.3 kpc and ra = 1, in four energy bands, compared 
to data (orange squares). The two cases ij = 14 and i? = 15 km give fully 
consistent results with ij = 13 km. Each profile is computed with the best 
fit viewing angles as reported in Table 2. 



not surprising, since the warm spot appears to be the remnant of 
the heated regions during the outburst. The elongated shape of the 
contour plots shows that the two angles ip and ^ are highly corre- 
lated in the fit. This is a result of the fact that the PF depends on a 
combination of these two angles. 

Having assumed a face-on spectrum to derive the spectral 
model parameters used to compute the pulsed fractions, we now 
proceeded to refit the spectrum, for each radius, using the best 
viewing geometry -0*,^*. This is iterated to a convergence crite- 



ria set by the measurement errors. The final spectral values in all 
cases are consistent with those found for i/i = ^ = 0. (Table 1). 
For example, for the special cases of R = 13 km, D = 3.3 kpc, 
and V* = 38° and C = 38°, we find hh = 7.5 ± 0.1, 
/cTbb = 0.453 ±0.003, fcTp = 0.218 ±0.001 and/?^ = 2.0±8:?, 
with x^ ~ 87.7 for 75 d.o.f. (see Figure|7j; as another example, for 
R=nkm,D = 2.0 kpc, and V* = 42° and C = 28°, we find 
riH = 8.4 ± 0.1, fcTBB = 0.413 ± 0.003, fcTp = 0.24 ± 0.03 and 
Ph = 2.0±?):?, with x^ = 76.4 (for 75 d.o.f.). Given the consis- 
tency between the spectral parameters determined with ^ = £, — 
and with those for the best viewing angles i/i* , ^* , the iterative pro- 
cess did not need to be continued further. For the best fit spectral 
parameters and viewing angles, the resulting pulsed fractions and 
pulse profile, for several radii, are show in Figure 3 and 5, respec- 
tively. It is iinportant to note that, for a pre-determined emission 
pattern of the radiation (here parameterized by the beaming param- 
eter n), there are values of radii (e.g., R ^ 12 km for n — 1) for 
which no good solution to the combined spectral/timing properties 
of the source can be found. Hence this type of analysis provides 
useful limits on the NS radius in the context of the model. 

To consider the effect of the angular distribution of radia- 
tion on our spectral and timing modeling, we repeated the above 
analysis for various assumed beaming indexes. Table |3]reports the 
model predicted PF for ^ n $: 2.5, in 0.5 step increments, for 
R = 9 km and R = 15 km (see Figure|6j. For n — (isotropic 
emission), we found that not even the largest radii were able to ac- 
count for the observed level of modulation. The modeled PF is, in 
fact, produced by the interplay of two flux components: the one 
from the whole NS surface, and the other from the small, hot spot. 
The high level of pulsation which could have been produced by the 
hot component alone, is strongly reduced by the presence of emis- 
sion on the entire NS surface. If the local emission is isotropic, not 
even the largest radii can provide sufficiently large PFs. However, 
as the value of n increases, the range of allowed radii becomes 
wider; we found that for n = 0.5, only R = 15 km provides a 
marginal match for the energy-dependent PF; for n = 1, the range 
R ^ 13 km provides a good overall representation of the data. For 
n = 2.5, due to the resulting strong increase of the PF at all ener- 
gies, a good match to the PFs can be found for each value of R (but 
clearly for different combinations of the angles i/;, ^). While our re- 
sults provide a hint to the presence of a magnetized atmosphere on 
the NS surface, they also show that without a priori knowledge of 
f(5) for each (6, <j)) on the NS surface, the timing analysis does not 
allow to break the degeneracy with the radius (at least within our 
measurement uncertainties). 



6 SUMMARY AND DISCUSSION 

The decline of the post-outburst emission of XTE J1810— 197 has 
allowed us to perform a detailed spectral and timing analysis of the 
emission of this source in its quiescent state. The main results from 
the modeling described herein can be summarized as follows: 

(i) The temperature distribution on the surface of the NS is con- 
sistent with the expectations of a dipole magnetic field configura- 
tion. We remark that our analysis does not guarantee that a dipo- 
lar magnetic field represents the only possible solution; in princi- 
ple, other magnetic field configurations and beaming patterns can 
be ad-hoc produced and tested for this source. However, a dipo- 
lar field is the simplest magnetic field configuration for a pulsar, 
and one that yields a single peak, nearly sinusoidal (as observed 
in XTE J1810— 197), for a wide range of viewing angles. Hence, 
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Table 2. Best emission geometry angles ^* and ^* and corresponding predicted pulsed fraction, for n 
uncertainty is reported 
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Table 3. Best emission geometry angles ^* and ^* and corresponding predicted pulsed fraction, for different value of the beaming factor n and D - 
compared to data. Ict c.l. uncertainty is reported 
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although not formally unique from a mathematical point of view, 
the proposed solution is physically motivated, and the fact that it 
provides such a good match to the data yields confidence that it is 
indeed a reasonably good representation of the quiescent emission 
of this transient magnetar. 

(ii) The pulsed fraction of the NS emission requires an 
anisotropic, pencil-type radiation pattern, which is a likely indica- 
tion of the presence of a magnetized atmosphere on the NS surface. 
For a cos" 5 emission profile, we performed a timing analysis for 



n in the range 0-2.5 (in steps of 0.5), and radii between 9-15 km 
(in steps of 1 km). We found that no match to the PF data could 
be obtained for n = (isotropic radiation pattern), no matter the 
value of the radius, while for n = 2.5, a good match to the PFs 
could be found even for R = 9 km (smaller values of the radius 
are increasingly allowed as n increases). For each value of n in be- 
tween there is a range of radii which are not allowed by the PF data 
(e.g. for n = 1 no good match could be found for _R ^ 12 km). 
Therefore, our analysis has clearly demonstrated how a detailed a 



The quiescent emission ofXTE J1810—197 L9 







* 


+ 
^ 








- 




ii 


40 


- 










_ . _ . _ 






20 


n=2.5 
n=2.0 

n=1.0 

n=0.5 
n=0.0 


























, 


II 





1 



40 



20 - 







1 

1 

t 1 




- 


- 


^~-.~ 


n=2.5 
n=2.0 
11=1.5 
n=1.0 
11=0.5 
n=0.0 









2.5 



1.5 2 2.5 1 1.5 2 

Energy [keV] Energy [keV] 

Figure 6. Model predicted PF (starting from the best emission geometry as reported in Table|3) for different values of the beaming pattern n. Left panel is for 
R=9 km, while right panel is for R=15 km. Icr c.l. uncertainty is reported. 



priori knowledge of the spectral distribution and emission pattern 
as a function of {6, (j>) on the entire surface of the star has the power 
to allow a radius constraint for the NS. 

(ill) The overall emission geometry is constrained by identifying 
likelihood regions in the 4^ ~ ^ parameter space. The most signifi- 
cant range according to our fits is consistent with the results of the 
earlier, post-outburst analysis of Pema & Gotthelf (2008), but the 
3(T allowed region is now reduced by a factor of 2.5. We note that 
the hot-spot axis used in the earlier work coincides with the dipole 
axis in the current study. 

A spectral and timing analysis of the post-outburst emission of 
XTE J1810- 197 has been performed also by Albano et al. (2010). 
They fitted the first 7 observations using a three-temperature model. 
The X-ray emission is produced in a globally twisted magneto- 
sphere, and parameterized by the twist angle, the electron velocity 
and the seed photon temperature. In the the last two sets, the flux is 
approaching the quiescent level, and only two temperatures are re- 
quired to fit the X-ray spectrum. For these observations, they found 
that the angle between the magnetic axis and the rotation axis is 
^ — 30.0°_20.o in the 6th observation, and ^ — 22.7°l2o'.o ^^ the 
7th, 6 months later, while the angle between the line of sight and 
the rotation axis was found to be -0 = 153.9°!tig;o ^^ 'he 6th and 
Tp = 145.8° _9 5 in the 7th observation. While the values of the 
angle ^ are roughly within the range of what we find, the angle ^ is 
inconsistent with ours. A direct comparison between our modeling 
and theirs is however not possible due to some fundamental differ- 
ences in the basic assumptions. Our quiescent emission is assumed 
to be thermal, and we leave the anisotropy level of the local radi- 
ation as a free parameter. On the other hand, the twisted magneto- 
sphere model has a well determined angular radiation pattern. The 
NS surface temperature is assumed to follow a dipolar-like pattern 
in our model, while in theirs it is assumed to be constant. Further- 
more, our modeling includes the general relativistic effect of light 
deflection (which heavily influences the pulsed fractions and hence 
the viewing geometry that we infer), while their modeling does not. 

Our study has allowed us to establish some important prop- 
erties of the quiescent emission of XTE J1810— 197 . The spec- 
tral analysis has shown evidence for the presence of a very small 
hot spot, only a couple of degrees in size. This emission compo- 
nent dominates the X-ray flux above about 1 keV. Interestingly, 



a small hot region is also characteristic of the quiescent emission 
of another (transient) AXP, CXOU J164710.2--455216 (Israel et 
al. 2007; Skinner et al. 2006). The most natural interpretation, 
in the case of XTE J1810-197 as wefl as in the case of CXOU 
J164710. 2— 455216 , is the association of this warmer region to a 
concentration of magnetic field lines, since the surface heat flux of 
a magnetar scales sensitively as B** (Thompson & Duncan 1996). 
It is then not surprising that this region coincides with the location 
of the outburst, as our analysis for XTE J1810— 197 has shown. 

Our constraints on the viewing geometry could in principle be 
combined with those obtained from radio measurements following 
the outburst, to further restrict the allowed parameter space. Camilo 
et al. (2007a) showed that the radio and X-ray locations were 
aligned, within some rather large uncertainty due to the small num- 
ber of counts of their Chandra data. Camilo et al. (2007b) used ra- 
dio polarimetry to fit for the angle angle between the magnetic field 
and the rotation axis (^ in our notation), and for the angle between 
the magnetic field axis and the line of sight (amin = "i/) — ^ in our 
notation), finding two possible configurations: one with ^ = 70° 
and Umin ~ 20° — 25°, and another with ^ = 4° and Umin = 4°. 
Polarization observations of XTE J1810— 197 were also used by 
Kramer et al. (2007) to constrain the viewing geometry of this pul- 
sar. They identified a main pulse and a mid pulse, produced in dif- 
ferent locations on the star. Their best fitted angles (with the same 
notation as above) were ^ — 44° and Umin ~ 39° for the main 
pulse, and ^ = 76° and Umin = 6° for the inter pulse. Given 
some discrepancies between the different angle estimates from ra- 
dio data, and the lack of a firm and accurate association between the 
radio and X-ray emitting regions, we refrain in the current work 
from using the radio angle determinations to restrict our allowed 
parameters range. However, we take the opportunity to emphasize 
the importance of simultaneous radio-X ray observations should 
another outburst occur (either in this or in another magnetar). 

Our analysis has allowed us to explore the consistency of the 
data with the expectations for the temperature distribution on the 
NS surface if the B field is dominated by a dipolar component. 
The angles ^ and tp that provide a good match to the observed 
pulsed fractions (Fig. 3 and Table 2) also ensure that the profile 
remains single-peaked in the lowest energy band, dominated by the 
surface/dipolar component. Combinations of the beaming parame- 
ter and the NS radius can be found for which our model provides an 
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excellent match to spectra, pulsed fractions and, more generally, to 
the full pulse profile. Within the context of magnetars, being able 
to infer the magnetic topology is important for a number of rea- 
sons. Firstly, measurements of the magnetic field strength from P 
and P make the implicit assumption of a dipolar field; substantial 
departures from this configuration will result in biased estimates. 
Second, the magnetic field topology in magnetars also plays an im- 
portant role in the resulting outbursting behaviour (Thompson & 
Duncan 1995, 1996), in that it regulates the intensity, frequency 
and location on the NS surface of the starquakes, and hence of the 
outbursts (Pema & Pons 2011). 

Ultimately, we note how the combined spectral/timing analy- 
sis that we performed has the potential to constrain the compact- 
ness ratio of the neutron star (or equivalently the radius, for a given 
mass). However, obtaining this type of constraint requires an a pri- 
ori knowledge of the spectral and radiation pattern of the local 
emission. This is particularly difficult in the case of highly magne- 
tized objects, since the local emission depends on the local strength 
and orientation of the B field, which in turn would need to be de- 
termined as a part of the fitting procedure. This analysis, which has 
begun for NSs with magnetic fields in the 10^^ — 10^^ G range 
(Ho & Mori 2008), will become possible for magnetars once atmo- 
spheric models for arbitrary field orientations (Lloyd 2003a,b; Ho 
et al. 2007) are extended to strengths B ~ lO" - 10^^ G. 
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Figure 7. Upper panel: XMM-Newton phase-averaged spectrum of XTE J 1810— 197 obtained with the model presented in the text, for the specific case of 
R=13 km, D=3.3 kpc, and the best fit angles ip* = 38° and ^* = 38°. Model residuals are shown in the lower panel. The parameters for the best fit model are 
reported in Table[T] The data show the summed spectra from the three observations. Lower panel: Same as the upper panel, but for unfolded source spectra. 
The contribution of the different model component is also shown: surface (black dotted line), warm spot (red dotted line). 



